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Probing the QCD vacuum with overlap fermions * 
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We use low lying eigenvectors of the overlap-Dirac operator as a probe of the QCD vacuum. If instantons 
play a significant role one would expect the low lying eigenmodes of the overlap-Dirac operator to consist mainly 
of the mixed "would be zero modes". Then, the eigenmodes should exhibit local chirality. Studying a recently 
introduced local chirality parameter, we find evidence supporting this picture. 



1. Introduction 

Understanding the QCD vacuum has been a 
goal of lattice QCD for many years. This under- 
standing, in turn, should illuminate the mecha- 
nisms of confinement and chiral symmetry break- 
ing. A good candidate for the source of chi- 
ral symmetry breaking are topological fluctu- 
ations of the gauge fields associated with in- 
stantons and anti-instantons. Through mixing, 
the "would be zero modes" from the instantons 
and anti-instantons become the small modes that 
lead, when present with a finite density, to chi- 
ral symmetry breaking, since {ipij}) = 7rp(0 + ), 
with /o(0 + ) the density of the small modes. Many 
of the low energy properties of QCD can be ex- 
plained phenomenologically by the interactions of 
the fermions with instantons and anti-instantons 
|IJ|^]. To our knowledge, confinement, though, is 
not one of those properties. 

There is a problem, though, with investigat- 
ing the topological content of Monte Carlo gen- 
erated gauge configurations. Ultraviolet fluctua- 
tions (noise) dominate these configurations. Pure 
gauge probes of the vacuum therefore need fil- 
ters, like smearing or cooling, to get rid of the 
UV noise and make it possible to uncover the 
underlying topological fluctuations ||. But do 
such filters change the physics? Fermions, on the 
other hand, should provide a more physical filter, 



provided lattice artifacts do not impede the chi- 
ral and topological properties. Since a few years 
such fermions are known, the overlap fermions of 
Neuberger [Q. 

Global topology, i.e. the topological suscepti- 
bility xt has already been successfully computed 
from the number and chirality of the zero modes. 
This was actually done, completely equivalently, 
from studying the spectral flow of the underlying 
Wilson Dirac operator 0. 

Here we want to address local properties. If in- 
stanton and anti-instanton like fluctuations are 
important, then the low- lying non-zero eigen- 
modes should have significant contributions from 
mixed "would-be zero modes" . Locally, in their 
peaks, these low eigenmodes should therefore be 
chiral, since zero modes are chiral, even though 
the integrated chirality of all non-zero modes is 
exactly zero. This can be measured with the "lo- 
cal chirality" parameter, X{x), defined by [p| 



tan 



(~(1 + X(x)) 



(1) 
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where ^Pl{x) and ipn(x) are the left- and right- 
handed components of the eigenmode. 

It should be noted, however, that Witten 
pointed out an inconsistency between instanton 
based phenomenology and large- N c QCD 0. In- 
stantons would produce an rjf mass that van- 
ishes exponentially for large N c , while consid- 
erations based on large- N c chiral dynamics sug- 
gest that the rf mass should be of order 1/N C . 
The topological charge fluctuations then should 
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be associated not to instantons but rather some 
other, confinement-related vacuum fluctuations. 
And the strong attraction produced by these 
confinement-related fluctuations would also in- 
duce the breaking of chiral symmetry. Then there 
is no reason that the low-lying modes should be 
chiral in their peaks. Using non-chiral Wilson 
fermions at a lattice spacing of a ~ 0.17 fm 
Horvath et al. || found that the local chirality 
X(x) was not peaked near ±1. We use chiral 
overlap fermions instead to avoid the potentially 
large 0(a) lattice artifacts from the explicit chiral 
symmetry breaking term of Wilson fermions. 

2. Setup 

For our study we consider the standard mass- 
less overlap-Dirac operator 

D(0) = ±[l + 75 e(H w (M))] (2) 

with H W (M) = jD w (-M) and D W (M) the 
usual Wilson-Dirac operator. Then, H 2 (0) = 
-D^(0)-D(0) commutes with 75 and can be simul- 
taneously diagonalized • H 2 (0) can have zero 
eigenvalues with chiral eigenmodes, which are 
also eigenmodes of D(0), due to global topology. 
The non-zero eigenmodes of H 2 (0) are doubly de- 
generate and have opposite chirality, 

# 2 (0)^U = , 75^T = >75VU = 

(3) 

with < A < 1. The non-zero (right) eigenmodes 
of D(0) are then easily obtained as 

1>± = -J=(l>l±*1>l) ( 4 ) 
with eigenvalues 

A± = A 2 ±iAVl-A 2 . (5) 

The "chirality" parameter of eq. ([!]) is then given 
by 

t <± [i+x{x)) ) = {ii^)) • (6) 



Table 1 

The gauge ensembles. Listed are the volume in 
lattice units and in units of fin , the number of 
configurations analyzed, the average number of 
zero modes and the average of the square of the 
global topological charge as obtained from the 
number of zero modes. 
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and is equal for the two related modes ip±. In 
the following they will be therefore counted as 
one mode. The exact zero modes, of course, have 
X(x) = ±1. 

We computed low- lying eigenmodes V^i of 

H 2 (Q) with the Ritz functional algorithm of 

Ref. |Hj. For the sign function e(H\y) in eq. (||) 
we used the optimal rational approximation of 
HH with "projection of low-lying eigenvectors of 
H w " to ensure sufficient accuracy. 

3. Results 

For a more extensive and detailed discussion of 
our results see plj |. Here, we concentrate on the 
results obtained on gauge field configurations gen- 
erated with the Iwasaki action |l^] . At compara- 
ble lattice spacing, the implementation of overlap 
fermions on such configurations is about a factor 3 
faster than on Wilson action configurations IO]. 
We computed the 20 eigenmodes i/'t,! 01 H^Jt)) 
with smallest A, with the Wilson-Dirac mass in 
the overlap-Dirac operator set to M = 1.65. Six 
ensembles were considered corresponding to three 
different lattice spacings and three different vol- 
umes, in physical units (see Table [l]). 

The chirality histograms of the two lowest non- 
zero modes at the 2.5% of sites with largest 
ip^ip(x) for all six ensembles are shown in Fig. [|. 
The physical volume of the systems shown in 
the left column are approximately the same, 2 
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Figure 1. Chirality histograms for the lowest two 
non-zero modes of the overlap-Dirac operator at 
the 2.5% sites with the largest ip^ip(x). 
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Figure 2. Chirality histogram for the lowest six 
non-zero modes at the 2.5% sites with the largest 
^4>{ x ) on the two lattices with volume of about 
7 fm 4 . 
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Figure 3. Chirality histogram for all the non- 
zero modes out of the 20 lowest modes at the 
20% sites with the largest ip^ip(x) on the lattice 
at /3 = 2.2872 with the largest volume. 



fm 4 . We see that the lowest non-zero modes 
tend to become more chiral in their peaks as the 
lattice spacing is decreased. This strongly sug- 
gests that the behavior observed here will sur- 
vive the continuum limit. Comparing the three 
histograms for (3 = 2.2872 having equal lattice 
spacing, we see a dramatic dependence on the 
physical volume. In an instanton liquid picture of 
the vacuum, the number of instantons and anti- 
instantons, and hence the number of almost zero 
modes grows linearly with the volume. Then it 
is not surprising that the lowest modes become 
increasingly chiral in their peaks with increasing 
volume. 

In Fig. H we show the chirality histogram of 
the lowest six non-zero modes on the two sys- 
tems with physical volume of about 7 fm 4 and 



in Fig. g the chirality histogram from all non- 
zero modes out of the 20 lowest modes that we 
computed (there are between 11 and 20 non-zero 
modes per configuration) for the system with the 
largest physical volume. Here we kept the 20% 
sites with largest %jpij){x). Comparing with the 
histograms in the left column of Fig. [j] indicates 
that the number of non-zero modes with similar 
chirality histograms grows roughly like the phys- 
ical volume so that there is a finite density of 
modes which are chiral in their peaks. This is 
good evidence that we are not observing a finite 
volume effect that would disappear in the infinite 
volume limit. 

In the chirality histograms shown we typically 
considered the contribution from the 2.5% sites 
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Table 2 

The average of the sum of tp'rp(x) over the 2.5% 
of sites kept for the chirality histograms as a per- 
centage of the sum over all sites for some low- lying 
eigenvectors for the various gauge field ensembles. 
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Figure 4. Chirality histogram for the lowest ten 
non-zero modes of the overlap-Dirac operator at 
the 6% sites with the largest ^ip(x) on 24 2 U(l) 
configurations at j3 = 1.9894. 



with the largest ip^ip(x). In Table || we list the 
sum of tp^tpfa) over those sites as a percentage of 
the sum over all sites for some low-lying eigen- 
vectors for all the gauge field ensembles consid- 
ered. The more chiral the eigenvectors are on 
those 2.5% sites, the higher is the percentage of 
the eigenvector contained on those sites. 

4. Abelian gauge theories 

To contrast the results for quenched QCD of 
the previous section we considered the chirality 
parameter also for quenched U(l) theories in two 
and four dimensions. In 2d topology plays a cru- 
cial role and one expects local chirality in the 
peaks of low-lying modes. This is clearly seen in 
Fig. ||, where the 10 lowest non-zero modes have 
been kept. 

In 4d U(l), in the confined phase, chiral sym- 
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Figure 5. Chirality histogram for the lowest two 
non-zero modes of the overlap-Dirac operator at 
the 1% (left panel) and 0.25% (right panel) sites 
with the largest tp^tp(x) on 8 4 U(l) configurations 
at = 0.9. 



metry is spontaneously broken, and a finite den- 
sity of near zero modes exists. But there are 
no instantons, and so one would not expect the 
near zero modes to be dramatically chiral even 
in their peak regions. There do exist exact zero 
modes of the overlap-Dirac operator, which again, 
of course, are chiral, although their origin is 
not completely understood fl3||. We analyzed, 
stored eigenmodes from Ref. [ p_3| . The chirality 
histogram for the lowest two non-zero modes is 
shown in Fig. |[ We note that while there is some 
mild indication of chirality peaking, the propor- 
tion of sites showing this behavior is much re- 
duced compared to the SU(3) case and appears 
to be of a qualitatively different nature than seen 
before. Also, the sum of ip^ip{x) over the sites 
kept for the histograms of Fig. [5] where almost 
identical for all 12 lowest eigenmodes (4.0% and 
1.20%, respectively; for the 2.5% of sites with 
largest ijj ^ijj(x) the sum is about 8.4% for all eigen- 
modes as compared to the non-abelian case in Ta- 
ble ||). In the U(l) case the (near) lack of peaking 
could conceivably come about from the scenario 
outlined by Horvath, et. al. - namely from the 
confinement inducing vacuum fluctuations. 

5. Conclusions 

By studying the local chirality parameter, in- 
troduced in Ref. m , for low-lying non-zero eigen- 
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modes of the overlap-Dirac operator, we investi- 
gated the topological content of the QCD vac- 
uum. Considering several volumes and lattice 
spacings, we found convincing evidence for chi- 
rality of the low lying modes in their peak re- 
gion. Our results give evidence that the number 
of modes which are chiral in their peaks grows 
linearly with the volume so that there is a finite 
density of such modes, and that the chirality of 
the modes becomes more pronounced as the lat- 
tice spacing is decreased. Therefore, our observa- 
tions should remain valid in the continuum limit. 
Our observations are consistent with the topolog- 
ical fluctuations being instanton dominated and 
the low-lying eigenmodes being the mixed "would 
be zero modes" due to the instantons and anti- 
instantons. 

Our findings confirm the results of Ref. Jl4| , 
where a different version of overlap fermions cou- 
pled to smoothened gauge fields were used, and 
extend them by a careful study of finite lattice 
spacing and volume effects. Our results, on the 
other hand, contradict the conclusions of Horvath 
et al. |6) who studied the chirality parameter with 
Wilson fermions at a large lattice spacing. That 
the results of Ref. ^| are strongly affected by lat- 
tice artifacts was confirmed in a recent paper of 
Hip et al. Jig] . These authors considered an im- 
proved chirality parameter for Wilson fermions, 
clover improved Wilson fermions, and lattices at 
smaller lattice spacing than Ref. |(|. With these 
reductions of lattice artifacts, they conclude that 
instanton dominance of topological charge fluc- 
tuations is not ruled out by the response of (im- 
proved) Wilson fermions. Together with the re- 
sults from overlap fermions presented here and 
the lack of significant chirality enhancement in 
the 4-d U(l) model where instantons should not 
exist, the case for instanton domination in 4-d 
SU(3) gauge theory, as measured by the chirality 
parameter, becomes even more compelling. How- 
ever, for an opposite point of view see Ref. jig ]. 
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